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Abstract 

We analyze the unforced and deterministically forced Burgers equation in 
the framework of the (diffusive) interpolating dynamics that solves the so- 
called Schrodinger boundary data problem for the random matter transport. 
This entails an exploration of the consistency conditions that allow to inter- 
pret dispersion of passive contaminants in the Burgers flow as a Markovian 
diffusion process. In general, the usage of a continuity equation dtp = —V{vp), 
where v = v{x, t) stands for the Burgers field and p is the density of trans- 
ported matter, is at variance with the explicit diffusion scenario. Under these 
circumstances, we give a complete characterisation of the diffusive transport 
that is governed by Burgers velocity fields. The result extends both to the ap- 
proximate description of the transport driven by an incompressible fluid and 
to motions in an infinitely compressible medium. Also, in conjunction with 
the Born statistical postulate in quantum theory, it pertains to the probabilis- 
tic (diffusive) counterpart of the Schrodinger picture quantum dynamics. 
We give a generalisation of this dynamical problem to cases governed by non- 
conservative force fields when it appears indispensable to relax the gradient 
velocity field assumption. The Hopf-Cole procedure has been appropriately 
generalised to yield solutions in that case. 
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1 Burgers velocity fields and the related stochas- 
tic transport processes 

The Burgers equation, [Q, ^, recently has acquired a considerable popularity in the 
variety of physical contexts, |@-||20|. An exhaustive discussion of its role in accoustic 
turbulence and gravitational contexts, where the emergence of shock pressure fronts 
is crucial, can be found in the monograph |17 . 



As is well known, the logarithmic Hopf-Cole transformation, 0, allows to replace 
the nonlinear problem ("nonlinear diffusion equation", [|l|) by a linear parabolic 
equation. Because of this equivalence all gradient-type solutions of the Burgers 
equation are known exactly. 

At the moment we shall preserve the gradient form restriction for Burgers velocity 
fields, but consider a more general form of the Burgers equation that accounts for 
an external force field F{x, t): 

dtv + iv-V)v = uAv + F{x,t) . (1) 

Let us mention that many of recent investigations were devoted to the analysis of 
curl V = solutions that are statistically relevant in view of the random initial data 
choice and/or inclusion of the random forcing term (the random potential in the 
related Parisi-Kardar equation, [pJ]]). 



However, irrespective of whether we need or need not the statistical input, an issue 
of matter transport driven by those nonlinear velocity fields needs the knowledge 
of an exact evolution of concentration and/or density fields. Much in the spirit of 
early hydrodynamical studies of advection and diffusion of passive tracers p2| , 
see also 

This particular issue is addressed in the present paper, under a simplifying assump- 
tion of nonrandom initial data and deterministic force fields. 

Following the traditional motivation (applicable both to incompressible and in- 
finitely compressible liquids, we regard the stochastic diffusion process as a 
primary phenomenon responsible for the emergence of (1) and thus justifying the 
"nonlinear diffusion equation" phrase in this context. 

Knowing the Burgers velocity fields, one is tempted to ask what is the particular 
dynamics (of matter or probability density fields) that is consistent with the chosen 
Burgers velocity field evolution. The corresponding passive scalar (tracer or con- 
taminant) advection- in-a- flow problem, [|I4|, [Tl|, |16|, is normally introduced through 
the parabolic dynamics: 

dtT + {v-\/)T = uAT (2) 
see e. g. also [^, |2^, |2^. For incompressible fluids, (2) coincides with the con- 



ventional Fokker-Planck equation for the diffusion process. This feature does not 
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persist in the compressible case. 

While looking for the stochastic implementation of the microscopic (molecular) dy- 
namics (2), ||2^, |TB|, it is assumed that the "diffusing scalar" (contaminant 



in the lore of early statistical turbulence models) obeys an Ito equation: 

dX{t) = v{x, t)dt + V2^dW{t) (3) 

X(0) = fo ^ X{t) = X 

where the given forced Burgers velocity field is perturbed by the noise term rep- 
resenting a molecular diffusion. In the, by now conventional, Ito representation of 
diffusion-type random variable X{t) one explicitly refers to the standard Brownian 
motion (e.g. the Wiener process) ^/2i'W{t), instead of the usually adopted formal 
white noise integral /q ff(s)(is, coming from the Langevin-type version of (3). 

Under these premises, while taking for granted that there is a diffusion process 
involved, we cannot view equations (l)-(3) as completely independent (disjoint) 
problems: the velocity field v cannot be quite arbitrarily inferred from (1) or any 
other velocity-defining equation without verifying the consistency conditions, which 
would allow to associate (2) and (3) with a well defined random dynamics, and 
Markovian diffusion in particular, |2^ . 

In connection with the usage of Burgers velocity fields (with or without external 
forcing) which in (3) clearly are intended to replace the standard forward drift of the 
would-be-involved Markov diffusion process, we have not found in the literature any 
attempt to resolve apparent contradictions arising if (2) and/or (3) are defined by 
means of (1). In particular, the usage of a continuity equation dtp = — V('yp), where 
V = v{x, t) stands for the Burgers field and p is the density of transported matter, 
is at variance with the explicit diffusion scenario. Also, an issue of the necessary 
correlation (cf. [|1^], Chap. 7. 3, devoted to the turbulent transport and the related 
dispersion of contaminants) between the probabilistic Fokker-Planck dynamics of 
the diffusing tracer, and this of the passive tracer (contaminant) concentration (2), 
has been left aside in the literature. 

Moreover, rather obvious hesitation could have been observed in attempts to 
establish the most appropriate matter transport rule, if (l)-(3) are adopted. De- 
pending on the particular phenomeno logical departure point, one either adopts the 
standard continuity equation, 0, that is certainly valid to a high degree of ac- 
curacy in the low viscosity limit (we refer to the standard terminology that comes 
from viscous fluid models; here, u stands for the diffusion constant) z/ J, of (l)-(3), 
but incorrect on mathematical grounds if there is a diffusion involved and simulta- 
neously a solution of (1) is interpreted as the respective current velocity of the flow: 
dtp{x,t) = —'V[v{x,t)p{x,t)] . Alternatively, following the white noise calculus tra- 
dition telling that the stochastic integral X(t) = Jqv{X{s), s)ds + jQff{s)ds implies 
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the Fokker-Planck equation, one adopts, [Q: dtp{x, t) = i/Ap(x, t) — V[v{x, t)p{x, t)] 
which is clearly problematic in view of the classic Mc Kean's discussion of the prop- 
agation of chaos for the Burgers equation, EBI, and the derivation of the 



stochastic "Burgers process" in this context: "the fun begins in trying to describe 
this Burgers motion as the path of a tagged molecule in an infinite bath of like 



molecules" , 27 



To put things on the solid ground, let us consider a Markovian diffusion process, 
which is characterised by the transition probability density (generally inhomoge- 
neous in space and time law of random displacements) p{y, s,x,t) , < s < t < T, 
and the probability density p{x,t) of its random variable X{t) , < t < T. The 
process is completely determined by these data. For clarity of discussion, we do not 
impose any spatial boundary restrictions, nor fix any concrete limiting value of T 
which, in principle, can be moved to infinity. 

The conditions valid for any e > 0: 

(a) there holds limtisjz^ I\y~x\>ePiy^ ■'^^ = 0^ 

(b) there exists a (forward) drift b{x, s) = Umtisjir^ I\y-x\<eiy ~ ^)p{x, s, y, t)d^y, 

(c) there exists a diffusion function (in our case it is simply a diffusion coefficient u) 
a{x, s) = limtis^J\g_s\<,{y - xfp{x, s, y, t)d^y, 



are conventionally interpreted to define a diffusion process, [^, g5| . Under suitable 
restrictions the function: 

g{x, s) = j p{x, s, y, T)g{y, T)(fy (4) 

satisfies the backward diffusion equation (notice that the minus sign appears, in 
comparison with (2)) 

-dsg{x, s) = uAg{x, s) + [b{x, s) ■ V]g{x, s) . (5) 

Let us point out that the validity of (5) is known to be a necessary condition for the 
existence of a Markov diffusion process, whose probability density p{x, t) is to obey 
the Fokker-Planck equation. Here, the new velocity field, named the forward drift 
of the process b{x,t), replaces the previously utilized Burgers field v{x,t))): 

dtp{x, t) = z/Ap(x, t) — V ■ [6(x, t)p(x, t)] . (6) 

The case of particular interest in the nonequilibrium statistical physics literature 
appears when p{y, s, x, t) is a fundamental solution of (5) with respect to variables 
y,s, ^ see however |^ for an alternative situation. Then, the transi- 
tion probability density satisfies also the Fokker-Planck equation in the remaining 
X, t pair of variables. Let us emphasize that these two equations form an adjoint 
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pair, referring to the slightly counterintuitive for physicists, although transparent 
for mathematicians, [^, 34, ^ 37 1, issue of time reversal of diffusion processes. 

After adjusting (3) to the present context, X{t) = /q b{X{s), s) ds + \^2i'W{t) we 
realise, ||38|, ^ that for any smooth function f{x,t) of the random variable 



X{t) the conditional expectation value: 



UmAtio^iJ Pix,t,y,t + At)f{y,t + At)d'y - f{x,t)] = {D+f){X{t),t) = (7) 

= (9i + (6-V) + z/A)/(f,t) 
X{t) = X 

— * — * 

determines the forward drift b{x, t) (if we set components of X instead of /) and 
allows to introduce the local field of (forward) accelerations associated with the 
diffusion process, which we constrain by demanding (see e.g. Refs. BR BBI, P? 
for prototypes of such dynamical constraints): 

{DlX){t) = {Dj){X{t),t) = {dtb+ {h- V)6 + v/\h){X{t),t) = F{X{t),t) (8) 



where, at the moment arbitrary, function F{x,t) may be interpreted as the exter- 
nal deterministic forcing applied to the diffusing system, In particular, if we 
assume that drifts remain gradient fields, curl 6 = 0, under the forcing, then those 
that are allowed by the prescribed choice of F{x,t) must fulfill the compatibility 
condition (notice the conspicuous absence of the standard Newtonian minus sign in 
this analogue of the second Newton law) 



F{x,t) = vn{x,t) 
n{x,t) = 2u[dt^ + i(|- + V-6)] . 



(9) 



This establishes the connection of the forward drift b{x,t) = 2z/V$(x, t) with the 
(Feynman-Kac, cf. [[ 
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potential ^l{x,t) of the chosen external force field. 
The latter connection, without invoking the Feynman-Kac formula, is frequently 
exploited in the theory of Smoluchowski-type diffusion processes, when the Fokker- 
Planck equation is transformed into the associated generalised diffusion equation. 

One of distinctive features of Markovian diffusion processes with the positive 
density p{x, t) is that the notion of the backward transition probability density 
p^{y, s,x,t) can be consistently introduced on each finite time interval, say < 
s <t<T: 

p{x, t)p^{y, s, f, t) = p{y, s, X, t)p{y, s) (10) 
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so that / p{y, s)p{y, s, x,t)d^y = p{x,t) and p{y,s) = J p^:{y, s,x,t)p{x,t)d^x. This 
allows to define (cf. ^ ^ for a discussion of these concepts in case of the 

most traditional Brownian motion and Smoluchowski-type diffusion processes) 

lim^tio^^[x- j pM-t- ^t,x,t)0y] = {D_X){t) = h,{X{t),t) (11) 

(D_/)(X(t),t) = (5, + {K. ■ V) - z/A)/(X(t),t) 

Accordingly, the backward version of the dynamical constraint imposed on the local 
acceleration field reads 

{Dim) = {Dim) = F{X{tlt) (12) 

where under the gradient-drift field assumption, curl 6^, = we have explicitly in- 
volved the forced Burgers equation (cf. (1)): 

dtk + {k ■ V)K - jyAK = F (13) 

Here, [^, in view of 6* = 6 — 2z/V/np, we deal with F{x,t) previously 

introduced in (9). A notable consequence is that the Fokker-Planck equation (6) 
can be transformed to an equivalent form of: 

dtp{x,t) = -uAp{x,t) -V[K{x,t)p{x,t)] (14) 

which however describes a density evolution in the reverse sense of time (!). 

At this point let us recall that equations (5) and (6) form a natural adjoint 
pair of equations that determine the Markovian diffusion process in the chosen time 
interval [0,T]. Clearly, an adjoint of (14), reads: 

dj{x, s) = uAfix, s) - [K{x, s) ■ V]/(f , s) (15) 

where: 

/(x, s) = J p^{y, 0, X, s)f{y, 0)d^y , (16) 

to be compared with (4), (5) and the previously mentioned passive scalar dynamics 
(2), see e.g. also p4|. Here, manifestly, the time evolution of the backward drift is 
governed by the Burgers equation, and the diffusion equation (15) is correlated (via 
the definition (10)) with the probability density evolution rule (14). 

This pair only can be consistently utilized if the diffusion proces is to be driven 
by forced (or unforced) Burgers velocity fields. Certainly, the continuity equation 
postulated to involve the Burgers field as the current velocity, does not hold true in 
this context. 
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Let us point out that the study of diffusion in the Burgers flow may begin from 
first solving the Burgers equation (12) for a chosen external force field, next specify- 
ing the probability density evolution (14), eventually ending with the corresponding 
"passive contaminant" concentration dynamics (15), (16). All that remains in per- 
fect agreement with the heuristic discussion of the concentration dynamics given in 
Ref. [jl6|. Chap. 7.3. where the "backward dispersion" problem with "time running 
backwards" was found necessary to predict the concentration. 

All that means that equations (l)-(3) can be reconciled in the framework set by 
(4)-(16). Then, the "nonlinear diffusion equation" does indeed refer to consistent 
stochastic diffusion processes. 

We are now at the point where the Burgers equation and the related matter 
transport can be consistently embedded in the general probabilistic framework of the 
so-called Schrodinger's boundary data (stochastic interpolation) problem, |31 



iD, n|, see also [^, In this setting, the familiar Hopf-Cole transformation. 



[5|, of the Burgers equation into the generalised diffusion equation (yielding 
explicit solutions in the unforced case) receives a useful generalisation. 



Indeed, in that framework, [g^, |3T|, the problem of deducing a suitable Markovian 
diffusion process was reduced to investigating the adjoint pairs of parabolic partial 
differential equations, like e.g. (5), (6) or (14), (15). In case of gradient drift fields 
this amounts to checking (this imposes limitations on the admissible force field 
potential, cf. also formula (9)) whether the Feynman-Kac kernel 

k{y, s, x,t) = J exp[- ^ c(w(r), r)dr](i/xj^''}(^) (17) 

is positive and continuous in the open space-time area of interest, and whether it 
gives rise to positive solutions of the adjoint pair of generalised heat equations: 

dtu{x, t) = iyAu{x, t) — c(x, t)u{x, t) (18) 

dtv{x, t) = —uAvi^x, t) + c{x, t)v{x, t) 

where c(x, t) = ^r2(x, t) follows from the previous formulas. In the above, (i/i||'jj(co') 
is the conditional Wiener measure over sample paths of the standard Brownian 
motion. 

Solutions of (18), upon suitable normalisation give rise to the Markovian diffusion 
process with the factorised probability density p{x,t) = u{x,t)v{x,t) wich interpo- 
lates between the boundary density data p{x,0) and p{x,T), with the forward and 
backward drifts of the process defined as follows: 

7>. X Vv(x,t) 

b{x,t)=2u—^-l- (19) 

V[X, t) 
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u{x, t) 

in the prescribed time interval [0,T]. The transition probabihty density of this 
process reads: 

p{y,s,x,t) = k{y,s,x,t) ,J . , (20) 

v{y,s) 

Here, neither k, (17), nor p, (20) need to be the fundamental solutions of appropriate 
parabolic equations, see e.g. ref. ||3l| where an issue of differentiability is analyzed. 

The corresponding (since p{x,t) is given) transition probability density, (10), of 
the backward process has the form: 

iiiy^ ts) 

P* {y, s, X, t) = k{y, s, X, t) J . (21) 

t) 

Obviously, |31|, |3^ , in the time interval < s < t < T there holds: 
m(x, t) = J uo{y)k{y, s, x, t)(Py and v{y, s) = J k{y, s, x, T)vt{x)(Px. 

By defining = logu, we immediately recover the traditional form of the 
Hopf-Cole transformation for Burgers velocity fields: 6* = — 2z/V$*. In the special 
case of the standard free Brownian motion, there holds b{x,t) = while h^{x,t) = 
—2uVlog p{x, t). 

Our discussion provides a complete identification of the stochastic diffusion pro- 
cess underlying both the deterministically forced Burgers velocity dynamics and the 
related matter transport, (14), the latter in terms of suitable density fields. The 
generalisation of the Hopf-Cole procedure to this case involves a powerful method- 
ology of Feynman-Kac kernel functions and yields exact formulas for solutions for 
the forced Burgers equation. Let us stress that the connection between the Burgers 
equation and the generalised (forward) heat equation is not merely a formal trick 
that generates solutions to the nonlinear problem. The forward equation (18) in 
fact carries a complete information about the implicit backward stochastic evolution 
, that is a Markov diffusion process for which the Burgers- velocity driven transport 
is appropriate. Notice that the transition probability density (21) obeys the familiar 
Chapman-Kolmogorov formula. If we wish to analyze a concrete density field gov- 
erned by this process, any two boundary density data p(x, 0) and p{x, T) allow to 
deduce the ultimate form of the (more traditional, forward) diffusion process (20), 
by means of the Schrodinger boundary data problem, Then, the adjoint 

pair of equations (18) gives all details of the dynamics, with (19)-(21) as a necessary 
consequence. On the other hand, the presented discussion implies a direct import 
of the shock-type matter density profiles to the general nonequilibrium statistical 
physics of diffusion-type processes. 
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2 The problem of nonconservative forcing of Burg- 
ers velocity fields 



By embedding the Burgers equation in the Schrodinger's interpolation framework, 
we could consistently handle random transport that is governed by gradient velocity 
fields and gradient-type external conservative forces. The natural question at this 
point is how to incorporate the the non-gradient (rotational for example) velocity 
fields and especially the nonconservative forces. This question may be addressed 
without reservations only in the context of the forced Burgers equation. Recall that 
the Hopf-Cole transformation is applicable only in case of gradient velocity fields. 
Moreover, the involved Schrodinger's interpolation framework extends the issue to 
the domain of nonequilibrium random phenomena, where standard Smoluchowski 
diffusions, [^, are normally discussed in the case of conservative force fields (and 
drifts in consequence). 



Remark: Strikingly, an investigation of typical nonconservative e.g. electromag- 
netically forced diffusions has not been much pursued in the literature, although an 
issue of deriving the Smoluchowski-Kramers equation (and possibly its large friction 
limit) from the Langevin-type equation for the charged Brownian particle in the gen- 
eral electromagnetic field has been relegated in Ref. ||50|| , 
of the innocent-looking exercise (sic !). On the other hand. 



Chap. 6.1 to the status 
the diffusion of realistic 
charges in dilute ionic solutions creates a number of additional difficulties due to the 
apparent Hall mobility in terms of mean currents induced by the electric field (once 
assumed to act upon the system), see e.g. pi] , p2| and [Q. In connection with the 
electromagnetic forcing of diffusing charges, the gradient field assumption imposes 
a severe limitation if we account for typical (nonzero circulation) features of the 
classical motion due to the Lorentz force, with or without the random perturbation 
component. The purely electric forcing is simpler to handle, since it has a definite 
gradient field realisation, see e.g. for a recent discussion of related issues. The 



major obstacle with respect to our previous (Section 1) discussion is that, if we wish 
to regard either the force F, (8), (12), or drifts b, 6* to have an electromagnetic 
origin, then necessarily we need to pass from conservative to non- conservative fields. 
This subject matter has not been significantly exploited so far in the nonequilibrium 
statistical physics literature. 



With this additional (to Burgersian per se) motivation, let us analyse how the 
gradient velocity field (and conservative force field) assumption can be relaxed and 
nonetheless the exact solutions to the Burgers equation can be obtained, both in 
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the unforced and forced cases, while involving the primoridal Markovian diffusion 
process scenario. 



It turns out the the crucial point of our previous discussion lies in a proper 
choice of the strictly positive and continuous, in an open space-time area, function 
k{y,s,x,t) which, if we wish to construct a Markov process, has to satisfy the 
Chapman-Kolmogorov (semigroup composition) equation. It has led us to consider 
a pair of adjoint partial differential equations, (18), as an alternative to either (5), 
(6) or (14), (15). 

The Feynman-Kac integration is predominantly utilised in the quantally ori- 
ented literature dealing with Schrodinger operators and their spectral properties. 
We shall exploit some of results of this well developed theory. The pertinent 
Feynman-Kac potential c{x,t) in (17), (18) is usually assumed to be a continuous 
and bounded from below function, but these restrictions can be substantially relaxed 
(unbounded functions are allowed in principle) if we wish to consider general Marko- 
vian diffusion processes and disregard an issue of the bound state spectrum and this 
of the ground state of the (selfadjoint) semigroup generator, |2^. Actually, what 
we need is merely that properties of c{x,t) allow for the kernel k, (17), that is pos- 
itive and continuous. This property is crucial for the Schrodinger boundary-data 
problem analysis. 

Taking for granted that suitable conditions are fulfilled, |^8|, ^ , we can immediately 
associate with equations (18) an integral kernel of the time-dependent semigroup 
(the exponential operator should be understood as time-ordered expression, since in 
general H{t) may not commute with H{t') for r ^ r'): 

k{y,s,x,t) = [exp{- H{T)dT)]{y,x) (23) 

where H{t) = —uA + c(r) is the pertinent semigroup generator. Then, by the 
Feynman-Kac formula, we get an expression (17) for the kernel, which in turn 



yields (19)-(22), see e.g. [^. As mentioned before, (20) combined with (17) sets a 
probabilistic connection between the Wiener measure corresponding to the standard 
Brownian motion with b{x, t) = 0, and this appropriate for the diffusion process with 
a nonvanishing drift b{x, t), curl 6 = 0. 

Our main purpose is to generalise (23), so that the positive and continuous (semi- 
group) kernel function can be associated with stochastic diffusion processes, whose 
drifts are no longer gradient-fields. In particular, the forcing is to be nonconservative. 

Since we have no particular hints towards Feynman-Kac type analysis of rota- 
tional motions, it seems instructive to invoke the framework of the Onsager-Machlup 
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approach towards an identification of most probable patlis associated witli tlie un- 



derlying diffusion process, |55, Bq, 57|. In tliis context, tlie non-conservative model 



system has been investigated in Ref. Namely, an effectively two-dimensional 



Brownian motion was analyzed, whose three-dimensional forward drift b{x), b-^ = 
in view of d^bi 7^ dyb2 has curl 6 7^ 0. Then, by the standard variational argument 



with respect to the Wiener-Onsager-Machlup action, |56|, ^ 

1 /"^^ 1 

I{L{S,x,ty,h,t2} = — {-[S~b{x,t)]^ + iyVb{x,t)}dt , (24) 

Zu Jti I 

the most probable trajectory, about which major contributions from (weighted) 
Brownian paths are concentrated, was found to be a solution of the Euler-Lagrange 
equations, which are formally identical to the equations of motion 

'q,, = E^ X B (25) 

of a classical particle of unit mass and unit charge moving in an electric field E and 
the magnetic field B. The electric field (to be compared with Eq. (9)) is given by: 

E = -V$ (26) 
$ = _i(62 + 2z/V6) 

while the magnetic field has the only nonvanishing component in the z-direction of 

B = curl b = {0, 0, d^b2 - dybi} , (27) 

Clearly, B = curl A where A=b is the electromagnetic vector potential. The simplest 
example is a notorious constant magnetic field defined by 61 (x) = — -fx2, b2{x) = 

One immediately realizes that the Fokker-Planck equation in this case is incom- 
patible with traditional intuitions underlying the Smoluchowski-drift identification: 
the forward drift is not proportional to an external force, but to an electromag- 
netic potential. Nevertheless, the variational information drawn from the Onsager- 
Machlup Lagrangian involves the Lorentz force-driven trajectory. Hence, some prin- 
cipal effects of the electromagnetic forcing are present in the diffusing system, whose 
drifts display an "unphysical" (gauge dependent) form. 

On the other hand, if we accept this "unphysical" random motion to yield the 
representation with the nongradient drift A: dX{t) = A(X(t),t)dt + ^/2i'dW{t), and 
consider the corresponding pair (5), (6) of adjoint diffusion equations with A{x,t) 
replacing b{x,t), then (8) tells us that 

{DlX){t) =dtA+ {A-V)A + uAA = -—{xi,X2,0} = -E{x) (28) 
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where E{x) = ^{xi, X2,0}, if calculated from (26). 

We thus arrive at the purely electric forcing with reversed sign (if compared with 
that coming from the Onsager-Machlup argument, (26)) and somewhat surprisingly, 
there is no impact of the previously discussed magnetic motion on the level of dy- 
namical constraints (8), (13). The adopted recipe is thus incapable of producing 
the magnetically forced diffusion process that conforms with arguments of Section 
1. Our toy model is inappropriate and more sophisticated route must be adopted. 



In below, we shall invoke the Feynman-Kac kernel idea (23), This approach 
has a clear advantage of elucidating the generic issues that hamper attempts to de- 
scribe the diffussion processes governed by nonconservative (and electromagnetic in 
particular) force fields. 

The Burgers equation and the problem of its nongradient solutions will appear resid- 
ually byproduct of the more general discussion. 



Usually, the selfadjoint semigroup generators attract the attention of physicists 
in connection with the Feynman-Kac formula. Since electromagnetic fields provide 
most conventional examples of nonconservative forces, we shall concentrate on their 
impact on random dynamics. 

A typical route towards incorporating electromagnetism comes from quantal moti- 
vations via the minimal electromagnetic coupling recipe which preserves the selfad- 
jointness of the generator (Hamiltonian of the system). As such, it constitutes a part 
of the general theory of Schrodinger operators. A rigorous study of operators of the 
form — A + V has become a well developed mathematical discipline, |Q. The study 
of Schrodinger operators with magnetic fields, typically of the form —(V — iA)'^ + V, 
is less advanced, although specialised chapters on the magnetic field issue can be 



found in monographs devoted to functional integration methods, |^8|, mostly in 
reference to seminal papers |60, |61[ . 



From the mathematical point of view, it is desirable to deal with magnetic fields 
that go to zero at infinity, which is certainly acceptable on physical grounds as well. 
The constant magnetic field (see e.g. our previous considerations) does not meet 
this requirement, and its notorious usage in the literature makes us (at the moment) 
to decline the asymptotic assumption and inevitably fall into a number of serious 
complications. 

One obvious obstacle can be seen immediately by taking advantage of the ex- 
isting results, [^. Namely, an explicit expression for the Feynman-Kac kernel in a 
constant magnetic field, introduced through the the minimal electromagnetic cou- 
pling assumption H{A) = — |(V — iAY, is available (up to irrelevant dimensional 
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constants): 

e.p[-</f(.4)l(f..T = |5-|p5l(^)"^ (29) 

xexp{-^{x3 - y-if - ^coth{^t)[{x2 - 1/2)^ + (a^i - yif] - i^{xiy2 - X2yi)} ■ 

Clearly, it is not real (hence non-positive and directly at variance with the major 
demand in the Schrodinger interpolation problem, as outlined in Section 1), except 
for directions y that are parallel to a chosen x. 

Consequently, a bulk of the well developed mathematical theory is of no use for 
our purposes and new techniques must be developed for a consistent description 
of the electromagnetically forced diffusion processes along the lines of Section 1, 
i.e. within the framework of Schrodinger's interpolation problem. Also, another 
approach is necessary to generate solutions of the Burgers equation that are not in 
the gradient form. 



3 Forcing via Feynman-Kac semigroups 

The conditional Wiener measure (i/i|?'^j(a;), appearing in the Feynman-Kac kernel 
definition (17), if unweighted (set c(u;(r),r) = 0) gives rise to the familiar heat 
kernel. This, in turn, induces the Wiener measure Pw of the set of all sample 
paths, which originate from y at time s and terminate (can be located) in the Borel 
set A ^ B? after time t — s: Pw[^ = /a I '^/^(f tj('^) ~ Ia ^I^ where, for simplicity 
of notation, the (y, t — s) labels are omitted and stands for the heat kernel. 

Having defined an Ito diffusion X{t) = Jq b{x, u)du + \/2i'W{t) we are interested 
in the analogous path measure: -Px[v4] = Jy^dx J dfi^(^l}^{u}-g) = Jy^dfi{X). 

Under suitable (stochastic, [^) integrability conditions imposed on the forward 
drift, we have granted the absolute continuity Px ^ Pw of measures, which implies 
the existence of a strictly positive Radon-Nikodym density. Its canonical Cameron- 



Mart in- Girsanov form, IHSl 321, reads: 



'^^^^\y,s,x,t) = exp^[fh{X{u),u)dX{u) - \ f [h{X{u),u)fdu] . (30) 



dfi 2u Js 2 



s 



If we assume that drifts are gradient fields, curl 6 = 0, then the Ito formula 
allows to reduce, otherwise troublesome, stochastic integration in the exponent of 
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(30), to ordinary Lebesgue integrals: 



^ ^'b{X{u),u)dX{u) = ^X{t),t)-^X{s),s)- Pdu[dt^ + lv-b]{X{u),u) . 

Zv J s J s Z 

(31) 

After inserting (31) to (30) and next integrating with respect to the conditional 
Wiener measure, on account of (9) we arrive at the standard form of the Feynman- 
Kac kernel (17). Notice that (31) establishes a probabilistic basis for logarithmic 
transformations (19) of forward and backward drifts: b = 2v'S/logv = 2i/V$, 
6* = —2v'S/logu = — 2i/V$*. The forward version is commonly used in connec- 
tion with the transformation of the Fokker-Planck equation into the generalised 



heat equation, P6| , p^ . The backward version is the Hopf-Cole transformation, 
mentioned in Section 1, used to map the Burgers equation into the very same gen- 
eralised heat equation as in the previous case, 0, 

However, presently we are interested in non-conservative drift fields, curlb ^ 0, 
and in that case the stochastic integral in (30) is the major source of computational 



difficulties, H3, E^, for nontrivial vector potential field configurations. It explains 



the virtual absence of magnetically forced diffusion problems in the nonequilibrium 
statistical physics literature. 

At this point, some steps of the analysis performed in Ref. |^5|] in the context 
of the "Euclidean quantum mechanics", cf. also are extremely useful. Let us 



emphasize that electromagnetic fields we utilize, are always meant to be ordinary 
Maxwell fields with no Euclidean connotations (see e.g. Chap. 9 of Ref. |^ for the 
Euclidean version of Maxwell theory) . 

Let us consider a gradient drift-field diffusion problem according to Section 1, 
with (17), (31) involved and thus an adjoint pair (18) of parabolic equations com- 
pletely defining the Markovian diffusion process. Furthermore, let A{x) be the 
time-independent vector potential for the Maxwellian magnetic field B = curl A. 
We pass from the gradient realisation of drifts to the new one, generalizing (19), 
for which the following decomposition into the gradient and nonconservative part is 
valid: 

b{x, t) = 2u Wlog $(x, t) - A{x) , (32) 

We denote 9{x,t)=exp[^{x,t)] and admit that (32) is a forward drift of an Ito 
diffusion process with a stochastic differential dX(t) = [2z/^ — A]dt + \/2h'dW{t). 
On purely formal grounds, we deal here with an example of the Cameron-Martin- 
Girsanov transformation of the forward drift of a given Markovian diffusion process 
and we are entitled to ask for a corresponding measure transformation, (30). 

To this end, let us furthermore assume that 9{x, t) = 6 solves a partial differential 
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equation 

dt9 = -z/[V - ^A{x)]'^e + c{x, t)e (33) 

with the notation c(x, t) = t) patterned after (9). Then, by using the Ito 

calculus and (32), (33) on the way, see e.g. Ref. we can rewrite (30) as follows: 

^(^>,x,t)= (34) 
djjL 

= e^^' ^[/(2^^ - A){X{u),u)dX{u) - \ /(2z/^ - A)\X{u),u) du\ 

= §frr4^^^[-^ f[A{u)dX{u) + uiy ■ A){X{u))du + n{X{u),u)du]] , 
0[X[s),s) zv Js 

where X{s) = y,X{t) = x. 

More significant observation is that the Radon- Nikodym density (34), if inte- 
grated with respect to the conditional Wiener measure, gives rise to the Feynman- 
Kac kernel (23) of the non- self adjoint semigroup (suitable integrability conditions 
need to be respected here as well, [§^]), with the generator = — z/[V— ^A{x)]'^ + 
c{x,t) defined by the right-hand-side of (33): 

dte{x,t) = H^e{x,t) = 

= [-Z/A + A{x) • V + i(V ■ A{x)) - ^[A{^)? + c{x,t)]9ix,t) (35) 
= -uAOix, t) + A{x) ■ Ve{x, t)+c xix, t)e{x, i) . 

Here: 

CA{x,t) = c{x,t) + i(Vi)(f) - j^[A{^)? . (36) 
An adjoint parabolic partner of (35) reads: 

dt9, = -H*^9, = uAe, + V[A{x)9,] - ca{x, t)d, = (37) 

= z/[V + ^A{x)fd, - c(f, t)e, . 

Consequently, our assumptions (32), (33) involve a generalization of the adjoint 
parabolic system (18) to a new adjoint one comprising (33), (37). Obviously, the 
original form of (18) is immediately restored by setting A = 0, and executing obvious 
replacements 9^ ^ u, 9 ^ v. 
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Let us emphasize again, that in contrast to Ref. [Q, where the non-Hermitean 
generator 2vH^^ (33), has been introduced as "the Euchdean version of the Hamil- 
tonian" H = — 2z/^(V — ^^)^ + ^, our electromagnetic fields stand for solutions of 
the usual Maxwell equations and are not Euclidean at all. 

As long as the coefficient functions (both additive and multiplicative) of the 
adjoint parabolic system (35), (37) are not specified, we remain within a general 
theory of positive solutions for parabolic equations with unbounded coefficients (of 
particular importance, if we do not impose any asymptotic fall off restrictions), 
pOl , |6^ , |6^, |6^. The fundamental solutions, if their existence can be granted, usually 
live on space-time strips, and generally do not admit unbounded time intervals. We 
shall disregard these issues at the moment, and assume the existence of fundamental 
solutions without any reservations. 

By exploiting the rules of functional (Malliavin, variational) calculus, under an 
assumption that we deal with a diffusion (in fact, Bernstein) process associated 
with an adjoint pair (35), (36), it has been shown in Ref. |6^ that if the forward 
conditonal derivatives of the process exist, then (D+X)(t) = - A = b{x,t), 

(32) and: 

{DlX){t) = {D+X){t) X curlA{x) + Vfi(x,t) + ucurl {curl A{x)) (38) 

where X(0) = 0, X{t) = x, x denotes the vector product in and 2z/c = Q. 

Since B = curl A = fioH, we identify in the above the standard Maxwell equation 
for curl H comprising magnetic effects of electric currents in the system: curl B = 

fio[D + cqE + Jext] where D = e^E while Jext represents external electric currents. In 
case of = 0, the external currents only would be relevant. A demand curl curl A = 
V(Vy4) — A A = corresponds to a total absence of such currents, and the Coulomb 
gauge choice VA = would leave us with harmonic functions A{x). 

Consequently, a correct expression for the magnetically implemented Lorentz 
force has appeared on the right-hand-side of the forward acceleration formula (38), 
with the forward drift (32) replacing the classical particle velocity q of the classical 
formula (25). 

The above discussion implicitly involves quite sophisticated mathematics, hence 
it is instructive to see that we can bypass the apparent complications by directly 
invoking the universal definitions (7) and (11) of conditional expectation values, 
that are based on exploitation of the Ito formula only. Obviously, under an assump- 
tion that the Markovian diffusion process with well defined transition probability 
densities p{y, s, x, t) and p*{y, s, x, t), does exist. 

We shall utilize an obvious generalization of canonical definitions (19) of both 
forward and backward drifts of the diffusion process defined by the adjoint parabolic 
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pair (18), as suggested by (32) with A = A{x): 

b = - a , K = -2z/^ - A . (39) 

We also demand that the corresponding adjoint equations (35), (37) are solved by 
9 and 6'^, respectively. 

Taking for granted that identities {D+X)(t) = b{x,t), X{t) = x and {D_X){t) = 
b^:{x,t) hold true, we can easily evaluate the forward and backward accelerations 
(substitute (39), and exploit the equations (35), (37)): 

{D+b){X{t),t) = dtb + {b ■ V)b + uAb = (40) 

= bx B + u curl B + VQ 

and 

{DX){X{t),t) = dtt + it ■ V)6* - uAb, = (41) 

= 6* X 5 - z/ curl S + V . 

Let us notice that the forward and backward acceleration formulas do not coin- 
cide as was the case before ( cf. Eq. (8), (12)). There is a definite time-asymmetry 
in the local description of the diffusion process in the presence of general magnetic 
fields, unless curl B = 0. The quantity which is explicitly time-reversal invariant 
can be easily introduced: 

v{x,t) = ^{b + t)ix,t)^ (42) 

^{Dl + D^_){X{t))=vxB + VQ . 

As yet there is no trace of Lorentzian electric forces, unless extracted from the term 
Vr2(x, t). This step we shall accomplish in Section 4. 

For a probability density 9^,9 = p of the related Markovian diffusion process, 
3^, we would have fulfilled both the Fokker-Planck and the continuity equations: 



dtp = uAp — V(6p) = — V('yp) = —uAp — V(6*p), as before (cf. Section 1). 

In the above, the equation (41) can be regarded as the Burgers equation with 
a general external magnetic (plus other external force contributions if necessary) 
forcing, and its definition is an outcome of the underlying mathematical structure 
related to the adjoint pair (33), (37) of parabolic equations. 

Our construction shows that solutions of the magnetically forced (hence nongra- 
dient) solution of the Burgers equation (41) are given in the form (39). In reverse. 
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the mere assumption about the decomposition of drifts, (39), into the gradient 
and nongradient part, imphes that the corresponding evolution equation (41) is the 
Burgers equation with the nonconservative forcing. The force term has a specific 
Lorentz form. Although we invoke electromagnetism, the decomposition (39) can be 
regarded to refer to an abstract nongradient component. In analogy to the previous 
Onsagcr-Machlup example, (24)- (28), the fictituous Lorentz force term would arise 
anyway. 



4 Schrodinger's interpolation in a constant mag- 
netic field and quantally inspired generalisations 

Presently, we shall confine our attention to the simplest case of a constant magnetic 
field, defined by the vector potential A — {— f X2, +f xi, 0}. Here, B — {0,0, S}, 
= and curl B = which significantly simplifies formulas (32)-(42). 

As emphasized before, most of our discussion was based on the existence assump- 
tion for fundamental solutions of the (adjoint) parabolic equations (33), (37). For 
magnetic fields, which do not vanish at spatial infinities (hence for our "simplest" 
choice), the situation becomes rather complicated. Namely, an expression for 

B^ 

c^{x, t) = c(f , t) - —{xl + xl) (43) 

includes a repulsive harmonic oscillator contribution. 

For the existence of a well defined Markovian diffusion process it appears neces- 
sary that a nonvanishing contribution from an unbounded from above c{x, t) would 
counterbalance the harmonic repulsion. 

To see that this must be the case, let us formally constrain 9{x,t) — exp{^{x,t)) to 
yield (in accordance with (9)) the identity: 

c(x, t) = dt^ + i/[V$]^ + i/A$ = . (44) 

Then, we deal with the simplest version of the adjoint system (35), (37) where, in 

— * 

view of V • A = = c, there holds: 

dtO = -i/[V - ^A]''9 = -uA9 + A-Ve - -^[AfO (45) 

dtO, = u[v + ^A]% = ^^d* + A-ve, + ^[A]% . 
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With our choice, curl A = {0, 0, B}, the equations (45) do not possess a funda- 
mental solution, which would be well defined for all {x,t) E x i?"*": everything 
because of the harmonic repulsion term in the forward parabolic equation. We can 
prove (this purely mathematical argument is not reproduced in the present paper) 
that the function k{y, s, x,t), : 

B 1 

k(y,s,x,t) = \ -T— ^y^^ (46) 

\ B B B 

xexp{- ^^^ _ {xz-y^f - —cot{—{t- s))[{x2-y2f + {xi-yif\ - —{xxy2-X2yi)} 

only when restricted to times t — s < ^ is an acceptable example of a unique positive 
(actually, positivity extends to times t — s< ^) fundamental solution of the system 
(44), (rescaled to yield i/ — > |). Here, formally, (46) can be obtained from the 

— * — * 

expression (29) by the replacement A — > —lA. 

An immediate insight into a harmonic repulsion obstacle can be achieved after 
an a; — |/-plane rotation of Cartesian coordinates: x[ = xicos^ut) — X2sin{ujt), x'^ = 
Xisin{ujt) + X2Cos{u!t), x'^ = x^, t' — t with uj — -q^. Then, equations (45) get 
transformed into an adjoint pair: 

dtS = -uA'9 - cu^x'^ + x'i)9 (47) 

dt'O, = uA'e, + u^x'^ + x'i)e, . 

Notice that the transformation a; — > ia; would replace repulsion in (47) by the 
harmonic attraction. On the other hand, we can get rid of the repulsive term by 
assuming that c{x, t), (43) does not identically vanish. For example, we can formally 
demand that, instead of (44), c(x,t) = +§;;ix'l + X2) plays the role of an electric 
potential. Then, harmonic attraction replaces repulsion in the final form of equations 
(35), (37). 

As a byproduct, we are given a transition probability density of the diffusion 
process governed by the adjoint system (cf. (28)): 

dte^-uAe + A-ve (48) 
dte,^iyAe, + A-ve, . 

with A — -^{—X2,xi,0}. Namely, by means of the previous x — y-plane rotation, 
(48) is transformed into a pair of time adjoint heat equations: 

dt'O = -uA'e , dfO^ = vA'e^ (49) 
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whose fundamental solution is the standard heat kernel. 

Finding explicit analytic solutions of rather involved equations (35), (37) is a 
formidable task on its own, in contrast to much simpler-unforced or conservatively 
forced dynamics issue. 

Interestingly, we can produce a number of examples by invoking the quantum 
Schrodinger dynamics. This quantum inspiration has been proved to be very useful 
in the past, ||36|, At this point, we shall follow the idea of Ref. |0 where 



the strategy developed for solving the Schrodinger boundary data problem has been 



applied to quantally induced stochastic processes (e.g. Nelson's diffusions, [38||). 
They were considered as a particular case of the general theory appropriate for 
nonequilibrium statistical physics processes as governed by the adjoint pair (18), 
and exclusively in conjunction with Born's statistical postulate in quantum theory. 

The Schrodinger picture quantum evolution is then consistently representable 
as a Markovian diffusion process. All that follows from the previously outlined 
Feynman-Kac kernel route, p5|, p3, ISO, pTl p3, HO, ET|, based on exploiting the 



adjoint pairs of parabolic equations. However, the respective semigroup theory has 
been developed for pure gradient drift fields, hence without reference to any impact 
of electromagnetism on the pertinent diffusion process. And electromagnetism is 
definitely ubiquitous in the world of quantum phenomena. 

Let us start from an ordinary Schrodinger equation for a charged particle in an 
arbitrary external electromagnetic field, in its standard dimensional form. To con- 
form with the previous notation let us absorb the charge e and mass m parameters 
in the definition of A{x) and the potential 0(x), e.g. we consider B instead of -^B 
and instead 0/m. Additionally, we set u instead of Then, we have: 

idtij{x, t) = -z/(V - ^AfHx, t) + ^(P{x)ij{x, t) . (50) 

The standard Madelung substitution ip = exp{R + iS) allows to introduce the 
real functions 9 = exp{R + S) and 6** = exp{R — S) instead of complex ones 
They are solutions of an adjoint parabolic system (35), (37), where the impact of 
(50) is encoded in a specific functional form of, otherwise arbitrary potential c{x,t): 

c{x,t) = ^^{x,t) = j-[2Q{x,t) - m] (51) 

Q(f,t) = 2u'^^j^^ = 2u'[ARix,t) + [VR{x,m . 

The quantum probability density p(x, t) = ip{x, t)ip(x, t) = 9(x, t)9^{x, t) displays 
a factorisation p = 99^^ in terms of solutions of adjoint parabolic equations, which we 
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recognize to be characteristic for probabilistic solutions (Markov diffusion processes) 



of the Schrodinger boundary data problem (cf. Section 1), |3T|, 0. It is easy 

to verify the validity of the Fokker-Planck equation whose forward drift has the form 
(38). Also, equations (40), (41) do follow with n = 2Q - (f). 

By defining E = — V0 (with utilised instead of ^0), we immediately arrive 
at the complete Lorentz force contribution in all acceleration formulas (before, we 
have used curl B = 0): 

dtb + (b ■ V)6 + uAb = bx B + E + vcurl B + 2VQ (52) 

dtK + {K ■ V)6* - uAK = K X B + E - vcurl B + 2VQ . 

Moreover, the velocity field named the current velocity of the flow, v = |(6 + 6^,), 
enters the familiar local conservation laws (see also ||32| for a discussion of how the 
"quantum potential" Q affects such laws in case of the standard Brownian motion 
and Smoluchowski-type diffusion processes) 



dp = -V{vp) (53) 

dtV + {v-V)v = vx B + E + VQ . 

A comparison with (33)-(43) shows that equations (50)-(53) can be regarded 
as the specialized version of the general external forcing problem with an explicit 
electromagnetic (Lorentz force-inducing) contribution and an arbitrary term of non- 
electromagnetic origin, which we denote by c{x, t) again. Obviously, c is represented 
in (51), by iQ(x,t). 

We have therefore arrived at the following ultimate generalization of the adjoint 
parabolic system (18), that encompasses the nonequilibrium statistical physics and 
essentially quantum evolutions on an equal footing (with no clear-cut discrimination 
between these options, as in Ref. |^T|) and gives rise to an external (Lorentz) 
electromagnetic forcing: 

dtOix, t) = [-z/(V - ^A)^ - i-0(f ) + c(f , t)]e{x, t) (54) 

9A(f,t) = + + ^(p{x) - c{x,t)]e,{x,t) . 

A subsequent generalisation encompassing time- dependent electromagnetic fields is 
immediate. 

The adjoint parabolic pair (54) of equations can thus be regarded to determine a 
Markovian diffusion process in exactly the same way as (18) did. If only a suitable 
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choice of vector and scalar potentials in (54) guarrantees a continuity and positivity 
of the involved semigroup kernel (take the Radon-Nikodym density of the form (34), 
with Q — s> —0 + Q , and integrate with respect to the conditional Wiener measure), 
then the mere knowledge of such integral kernel suffices for the implementation of 
steps (18)-(22), with u ^ 6^,, v ^ 6. To this end it is not at all necessary that 
k{x,s,y,t) is a fundamental solution of (54). A sufficient condition is that the 
semigroup kernel is a continuous (and positive) function. The kernel may not even 
be differentiable, see e.g. Ref. for a discussion of that issue which is typical for 
quantal situations. 

After adopting (54) as the principal dynamical ingredient of the electromagnet- 
ically forced Schrodinger interpolation, we must slightly adjust the emerging accel- 
eration formulas. Namely, they have the form (52) , but we need to replace 2Q{x, t) 
by, from now on arbitrary, potential Q{x,t) = 2h'c{x,t). The second equation in 
(53) also takes a new form: 

dtV + {v-V)v = vx B + E + Vin-Q) (55) 

see e.g. Ref. |^2| for more detailed explanation of this step. The presence in (54) of 
the density-dependent — VQ term finds its origin in the identity 6 — 6=,, = 2z/Vp(x, t) 
and is a necessary consequence of the involved (forced in the present case) Brownian 
motion, see e.g. [§6], 0, 



Finally, the second equation (52) with replacing 2Q is the most general form of 
the Burgers equation with an external forcing, where the electromagnetic (Lorentz 
force) contribution has been extracted for convenience. Solutions of this equa- 
tion must be sought for in the form (39), which generalizes the logarithmic Hopf- 
Cole transformation to non-gradient drift fields. Equations (54) are the associated 
parabolic partial differential (generalised heat) equations, that completely determine 
probabilistic solutions (Markovian diffusion processes) of the Schrodinger boundary 
data (interpolation) problem. In turn, for this particular random transport, the 
forced Burgers velocity fields play the role of backward drifts of the process. 



5 Outlook 

Our discussion, albeit motivated by the issue of diffusive matter transport that is 
consistently driven by Burgers velocity fields (this extends both to the compressible 
and incompressible cases), has little to do with classical fluids. The emergence of 
shock pressure fronts is more natural in the compressible situation. This shock profile 
possibility (inherent to the Burgers equation) has been imported to the nonequilib- 
rium statistical physics of random phenomena by exploring the idea of Schrodinger's 
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interpolation problem and revealing its connection with the Burgers dynamics. That 
has been the subject of Section 1. 

The next important result (a preliminary discussion of rotational Burgers fields can 



be found in Ref. |23[) amounts to relaxing the gradient-field assumption (that is 
crucial for the validity of the Hopf-Cole transformation). In Sections 2 and 3 we 
have analyzed the ways to generalise the Feynman-Kac kernel strategy so that the 
involved (drifts) velocity fields admit the nongradient form. Our analysis was per- 
fomed with rather explicit electromagnetic connotations. Equations (35) and (37) 
generalise the adjoint pair (18) to diffusion processes with nongradient drifts, (39). 
As follows from (41), the very presence of the nongradient term in the decompo- 
sition (39) implies that the corresponding evolution equation for the velocity field 
(backward drift of the process) is the Burgers equation with the specific Lorentz- 
type forcing. 

Section 4 extends the discussion to quantally implemented diffusion processes, where 
the minimal electromagnetic coupling is a celebrated recipe. This quantal motivation 
allows to arrive at the adjoint system (54), that incorporates an electric contribution 
and allows to define and solve the Burgers equation with the combined conservative 
and nonconservative (electromagnetic, in particular) forcing. Let us emphasise again 
that a transformation of the Burgers equation (whatever the force term is) into a 
generalised diffusion equation is not merely a formal linearisation trick. This, [|l|, 
"nonlinear diffusion equation" does indeed refer to a well defined stochastic diffusion 
process, but a complete information about its features is encoded in the involved 
parabolic equations. 
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